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ABSTRACT 

We present far-infrared maps and spectroscopy for a number of deeply embedded protostellar 
objects (CepE, HH211-MM, IC1396W, L1157, L1211, and RN015FIR) from data that 
we acquired with the ISO instruments PHOT and LWS. Several previously undetected deeply 
embedded sources are found in the vicinity of our targets. We determine temperatures and 
luminosities of seven objects and locate them on a l^boi-'^boi diagram - the equivalent to a 
Hertzsprung-Russell diagram for protostars. Their masses and ages, according to their loca- 
tion on tracks taken from our evolutionary model, are derived. L 1211 and CepE appear to 
be intermediate mass objects which will reach final masses of about 3Mq, while the other 
sources are in or below the solar mass range. The derived ages of 15000 to 30000 yr are 
consistent with their current Class state. A comparison of the luminosity of the associated 
outflows in the 1 -0S(1) line of molecular hydrogen with the source properties (bolometric 
luminosity, bolometric temperature, and envelope mass) of 16 Class sources shows no sta- 
tistically significant correlations. Nevertheless, the data are consistent with a scheme in which 
the outflow strength and protostar evolve simultaneously. We show that the relationship is par- 
tially disguised, however, by the local properties of the surrounding material, the extinction, 
and short-term flux variability. 
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1 INTRODUCTION 

In the earliest stages of star formation - the so-called Class 
and Class 1 phases - protostars are still deeply embedded in their 
parental molecular cloud cores. This material absorbs almost all of 
the emitted radiation of the star in the optical and the near-infrared. 
The spectral energy distribution (SED) of Class protostars peaks 
at about 100-160 pim, the wavelength of the maximum of a modi- 
fied blackbody at 30-80 K. Hence, direct observations of protostars 
have to be carried out in the far-infrared and in the (sub-)mm wave- 
length range. Sub-mm and millimeter observations of some of the 
source s investigated here have b een obtained e.g. by Lefloch et al. 

1997|) . Chini et al. <200ll) . Gi jeth e t al. 
and Gueth & Guilloteau jl999l) . 
1996 



199€ ), Ladd & Hodapp 
1997 ), Motte 8l Andre JioOll) 



The ISO satellite (Kessler et al. <1996l) ) with its PHOT in- 
strument had the capacity to measure the broad-band continuum in 
the far-infrared. Such observations, covering the peak region of the 
SED of protostars, help to yield some of the major properties of 
these objects such as their temperature, the sub-mm slope of their 



* Based on observations with ISO, an ESA project with instruments funded 
by ESA Member States (especially the PI countries: France, Germany, the 
Netherlands and the United Kingdom) and with the participation of ISAS 
and NASA. 



SED, the optical depth and the solid angle under which they emit. 
The latter two cannot be disentangled due to the limited spatial res- 
olution of the ISOPHOT instrument. With higher resolution obser- 
vations (e.g. SCUBA), however, we can independently determine 
the solid angle under which an object is seen and that way infer its 
optical depth. These parameters, together with the distance, enable 
us to calculate the total (Lboi) and sub-mm (Lsmm) luminosities 
of each object. We may then decide whether an object really is of 
Class or not by determining the hsmmfLtoi ratio (Andre et al. 
i2OO0t) ). Finally, by placing the inferred values on a temperature - 
bolometric luminosity diagram - the equival ent to a Hertzsprung- 
Russell diagram for protostars (Myers et al. il998l) ) - we are able 
for the first time to estimate the (model dependent) ages and masses 
of these sources directly. 

Bipolar outflows invariably accompany Class sources: 
strong inflow and outflow of material are concurrent. We thus wish 
to probe how the mass outflow rate is related to the mass accretion 
rate onto the protostar. The outflowing material interacts with the 
ambient medium through radiative shocks. Thus, the luminosity of 
the outflow may be correlated with some of the source properties 
(e.g. the bolometric source luminosity), which depend on the mass 
accretion rate. Therefore, we measured the luminosities of the out- 
flows of 16 ClassO sources in the 1 -OS(l) line of molecular hy- 
drogen. This is usually the strongest and easiest line to observe in 
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Table 1. Log of our ISOPHOT and LWS observations. 



Observation 
number 


Object 


a (12000) 


S (J2000) 


AOT 


teip [s] 


65903003 


RNO 15 FIR 


03 27 39 


+30 13 00 


PHT22 


670 


65903004 


RNO 15 FIR 


03 27 39 


+30 13 00 


PHT22 


620 


65903101 


HH211 


03 43 57 


+32 00 49 


PHT22 


670 


65903102 


HH211 


03 43 57 


+32 00 49 


PHT22 


620 


65201107 


HH211 


03 43 57 


+32 00 52 


LWSOI 


2268 


65902801 


HH 211 West 


03 43 57 


+32 01 04 


LWSOl 


3350 


66600502 


HH211 East 


03 43 59 


+32 00 36 


LWSOI 


2912 


46601429 


LI 157 


20 39 06 


+68 02 13 


LWSOl 


3390 


28200120 


L1157 


20 39 06 


+68 02 14 


LWSOI 


1958 


52902105 


LI157 


20 39 06 


+68 02 14 


PHT22 


668 


52902106 


LI 157 


20 39 06 


+68 02 14 


PHT22 


620 


54301407 


IC 1396W 


21 26 06 


+57 56 17 


PHT22 


668 


54301408 


IC 1396W 


21 26 06 


+57 56 17 


PHT22 


620 


56300709 


L 1211 


22 47 17 


+62 01 58 


PHT22 


670 


56300710 


L 1211 


22 47 17 


+62 01 58 


PHT22 


620 


56600912 


Cep E South 


23 03 13 


+61 41 56 


LWSOl 


1888 


56402111 


CepE 


23 03 13 


+61 42 27 


PHT22 


670 


56402112 


CepE 


23 03 13 


+61 42 27 


PHT22 


620 


56601113 


Cep E North 


23 03 13 


+61 42 59 


LWSOl 


1890 



near-infrared spectra of shocked molecular hydrogen, and due to 
the short cooling time of H2 it is a good tracer of the present in- 
teraction of the outflow with the surrounding material. These H2 
luminosities are then compared with various source properties to 
investigate possible correlations. 

Modelling of Class protostars remains in its infancy. 
Schemes now exist which yield evolutionary tracks, b ased o n re- 
lating gas accretion to the dusty envelope ( Myers et al. Jl998l) ) and 
jet thr ust to ga s accretion ( Bontemps et al. jl99a). Saraceno et al. 
\l99dt) . Smith Jl998l:l2000h . and Andre et al. l200(t\ ). We combine 
these schemes here in order to test if the simplest assumptions, such 
as a spherical envelope and a single accreting object, are feasible. 

In this paper, we first present our far-infrared ISO maps and 
spectroscopy, and then summarize the data analysis and how we de- 
rive temperatures and luminosities (Sect.|2j. In Sect.|3| we present 
our results, and comment on individual objects. A discussion of age 
and mass determination, and the general relationship to the out- 
flows is contained in Sect.|4| A framework within which the data 
can be interpreted is then put forward (Sect.|5j. 



2 OBSERVATIONS AND DATA ANALYSIS 

We used the ISO satellite to obtain ISOPHOT minimaps of six 
Class sources and LWS full grating spectra for three of them. All 
observations are listed in TableQ 

2.1 ISOPHOT data 

Minimaps were taken for six objects (CepE, HH21I-MM, 
IC 1396 W, L 1 157, L 121 1, and RNO 15 FIR) with ISOPHOT in 
its PHT22 mode by single pointing and moving of the telescope 
by one (ClOO) or half (C200) of a detector pixel. We used four 
filters (60, 100, 160, and 200 /xm). For 60 and 100 /xm, the ClOO 
detector (3x3 array of Ge:Ga) was used to create a 5 x 3 pixel min- 
imap with a pixel size of 45" x 46". The maps thus cover a field 
of view of 230"x 135". For the two longer wavelengths 7x3 mo- 
saics with a pixel size of 45" x 90" were obtained using the C200 
detector (2x2 array of stressed Ge:Ga), covering thus a field of 
view of 315" X 270". For details on the instrument and the used 
Astronomical Observing Templates (AOT) see the ISO Handbook, 



Table 2. Observation log of the NIR observations. The used telescopes, de- 
tectors and filters are listed. H2 indicates the narrow band filter, centred at 
the 1 - 0S(1) line of H2. The narrow band filter at a wavelength of 2.140 fim 
(continuum) is labeled with 2140. The number of images is separately indi- 
cated for each filter. In some cases the investigated objects fill only a part of 
the whole obtained mosaic (esp. HH212). The observing time is given per 
single image. 



Observatory 


Telescope 


Object 


Filter 


Number of 




Dale 


Detector 








images 




La Silla 


ESO/MP12.2-m 


HH24 


H2 


, K' 


12,3 


20,2 


Apr93 


IRAC2 


Ser-FIRSI 




23,4 


20,2 






VLA 1623 


Ha.K- 


24,4 


20, 2 


Calar Alto 


2.2 m 


LI448 


H2 


, K' 


32,33 


60, 3 


Jan94 


MAGIC 












Calar Alto 


2.2 m 


DR21 


H2 


. K' 


437, 78 


25, 25 


Sep94 


MAGIC 


L 1157 


H2 


. K' 


140, 13 


25, 100 


Calar Alto 


3.5 m 


CepE 


H2 


,2140 


52, 120 


30, 30 


Nov95 


MAGIC 


HH211-MM 


H2 


,2140 


60,59 


30, 30 






L 1157 


H2 


,2140 


44,71 


30, 30 


Manna Kea 


UH2.2-m 


L 1211 


H2 


, K' 


28,27 


200, 60 


Ang97 


QUIRC 












Calar Alto 


3.5 m 


Cep A 


H2 


,2140 


202, 368 


20,20 


Sep97 


MAGIC 


HH21I 


H2,2140 


32, 27 


20, 20 






L 1448 


H2 


, K' 


247, 253 


20,3 


Calar Alto 


1.2 m 


Cep A 


H2.K' 


64,32 


60, 15 


Nov98 


MAGIC 


CepE 


H2 


, K' 


34, 14 


60, 15 






HH212 


H2 


, K' 


3112, 1487 


60, 15 






L 1448 


H2 


, K' 


400, 149 


60, 15 






RNO 15 FIR 


H2 


, K' 


120, 60 


60, 15 


Calar Alto 


3.5 m 


CepE 


H2 




24 


30 


DecOO 


OMEGA PRIME 


HH211 


H2 




41 


30 






LI 157 


H2 




17 


30 






LI448 


H2 




42 


30 



Vol ume V : PHT — The Imaging Photo Polarimeter and Lemke et 
al. h99d) . The data were reduced with the ISOPHOT Interactive 
Analysis (PIA V9.1) software. 

Flux measurements in the ISOPHOT maps were carried out 
in two different ways: I) Point spread function (PSF) photometry 
using PSF fractions provided by Laureijs il999t) was done for the 
ClOO maps. We do not provide PSF photometry for the C200 de- 
tector since the given PSF fractions by Laureijs 1 19991) are only for 
the whole C200 pixel and our maps have a sampling of half a pixel 
in one direction. 2) "Aperture" photometry was obtained for all 
filters of both ClOO and C200 detectors. Here we attributed each 
pixel in the maps either to 'object' or to 'background' manually, 
then summed up both and subtracted 'background' from 'object' to 
obtain its flux. Since at 60 and 100 ^m, i.e. for the ClOO data, we 
were able to do photometry with both methods, we have a means of 
estimating the consistency of both. All measured fluxes, including 
the available IRAS fluxes of our objects, and the background level 
in the maps are provided in Table|3| 

2.2 LWS data 

For three objects (L 1 157, Cep E, and HH 21 1) we have full grating 
medium-resolution LWSOl scans, which cover a wavelength range 
from 43 to 196.9 /im with a resolving power between 150 and 300. 
See the ISO Handbook, Volume IV: LWS — The Long Wavelength 

^ http://www.iso.vilspa.esa.eS/manuals/HANDBOOK/V/pht_hb/ 
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Spectrometer^ and Clegg et al. il99d) for instruments and AOT 
details. We reduced the LWS data using standard pipeline 7. For 
deglitching and flux calibration and defringing of the spectra we 
employed the ISO Spectral Analyses Package (ISAP 1.6a). 



J (l -e"^) B{\,T)Tf{X)dX 
x=o 

oo 

/ Tf{\)d\ 



(4) 



2.3 Near-infrared H2 observations 

For the measurement of the luminosities of the outflows in the 1 - 
0S(1) line of molecular hydrogen at 2.122/im near-infrared im- 
ages were taken in several observing campaigns and at various tele- 
scopes. The complete list of all observations is provided in Table|2| 
We observed the objects in two filters to distinguish between line 
and continuum emission. Due to the angular size of the objects, the 
single images had to be arranged into large mosaics. All observing 
campaigns were (re)-reduced for consistency using own software 
based on the IRAF package DIMSUM. The whole procedure in- 
cludes flatfielding, cosmic ray hit removal and sky subtraction as 
well as re-centering and mosaicing. For a higher astrometric accu- 
racy we used all available stars in the field for the re-centering. 
The photometric calibration was achieved by the observation of 
faint near-infrared standards with an accuracy of 10%. For the flux 
measurements we subtracted the scaled continuum image from the 
emission line image to measure only the flux in the 1 - S(l) line 
of H2. 

Our images are being prepared for publication or are already 
published. Since in this paper we will only use the integrated H2 
line luminosities we do not reproduce the images here. The objects 
are dis cussed in the follo wing pa pers: RNO 15 FIR in Davis et al. 

1997t) and Rengel et al. ilOO'j) : HH211-MM in Eisloffel et al. 

2OO3I) : VLA 1623, L 1157 in D avis and Eisloffel 1199^) : L 121 1 
in Froebric h and E isloffel t004) : CepE in E isloffe l et al. jl996h . 
Smith et al. ^2003j); L 1448 region in Eisloffel |2OO0|) and Froebrich 
et al. ( 2002); HH 212, HH 24 in Froebrich et al. 62001.) and Eisloffel 
et al. j|2004l) : Ser-FIRSl in Eisloffel and Froebrich j20ol . 



2.4 Fit of the spectral energy distributions 

The observed broad-band continuum fluxes of our sources allow us 
to fit a SED to the measurements and to infer source properties (e.g. 
Tboi and Lboi)- To fit the SED we used Eq.Qfor the flux density S 
of our objects. 



S[Jy]/T,D. = (l-e"") ■ B{X,T) 



(1) 



B{X,T) is the Planck function, ES7 the solid angle of the source 
and T the optical depth, r is set as 



T = TIOO 



A 



100 fim 



(2) 



A is in jj,m, the optical depth at 100 ^m (rioo) is a free parameter, 
and 13 is the sub-mm slope of the SED. The lowest rms of the fit is 
obtained when the solid angle of the object is determined by 



sn = 



^ (A 5^)2 




(A 57)2 



(3) 



/ indicates the various used filters, 5*™ the flux measurements in 
these filters, and A the error of the measurements. Sj/ESl is 
calculated by 



' |http://www.iso.vilspa.esa.es/manuals/HANDBOOK/IV/lws-hb/| 



for each filter / separately using the filter transmission curves 
TfiX). 

To fit a graybody to the measured SEDs of each object, a grid 
of graybodys was computed (see Eq.0 in which we varied the three 
parameters T, tioo, and (3. We varied the temperature between 15 
and 80 K, in steps of 0.25 K, the optical depth at 100 fim from 0.09 
to 40, in logarithmic intervals of 1.5, and the sub-mm slope from 
0.0 to 3.0, in steps of 0.1. These graybodys were convolved with the 
filter curves of the used filter bands (see Eq.|4j. Then the solid angle 
ESI was determined by computing the deviation of the model points 
from the measurements and minimising this value (see Eq.|3}. Fi- 
nally the rms of the fit to the measurements was calculated (see 
Eq.|5l n indicates the number of filters) and the parameters leading 
to the minimal rms were selected. 



'1 V- {SfiT,l3,Two,m-Sfr 



f 



(ASf 



(5) 



We find that rioo has almost no influence on the shape of the 
graybody curve, but only on the absolute flux level, which on the 
other hand mainly depends on ESI. Thus, the values presented in 
Table|4]are determined by fixing rioo to unity. This restriction has 
no influence in the deduced parameters Ttoi and htoi , but the given 
source size Ef2 has no physical meaning. In Sect.|3|we present also 
graybody fits with rioo as a free parameter, in case that this im- 
proves the fit significantly. If sub-mm or millimetre observations 
yield source sizes, the optical depth at 100 ^m can be constrained. 

Another way to determine physical meaningful radii for the 
sources is to follow the assumptions of Myers et al. il998h . They 
adopt an optically thick envelope, a single power law dependence 
with the frequency of the emissivity and an envelo pe den sity pro- 
portional to Taking Figure 2 in Myers et al. i 19981) . we can 



determine the optical depth using log rioo ~ 



+ 4.28. 



We used the bolometric temperatures from Table|4|to derive rioo, 
and repeated the fit of the SED with this optical depth. This leads to 
a new solid angle of the source which is the size of the protostellar 
envelope where rioo has the correct value, according to the assump- 
tions of Myers et al. il998h . In first approximation the solid angle 
and the optical depth are connected by ESl^^ oc (1 — e^^^""). 
Thus, we can determine the radius of the envelope where rioo is 
unity (presented as Rioo in Table|4j. This radius is different from 
the envelope sizes obtained of optically thin emission by sub-mm 
or millimeter measurements (e.g. Motte & Andre |2001) and Chini 
et al. i200ll) ). 

The above described method to fit the SED was applied to 
both, the measurements at our four ISOPHOT wavelengths and all 
available data from Table|3|(except the IRAS 12 and 25 /im points; 
see below). The inferred object properties from the latter are listed 
in Table|4| We did not find significant changes in the fit param- 
eters between the two methods, except for L 1211 and partly for 
RNO 15 FIR. These differences are discussed in Sect. l3.5l and l3.6l 
respectively. All obtained results are discussed for each object sep- 
arately in Sect.|3| 

IRAS 12 and 25 nm points are excluded from the fit of the 
SED since recent works show, that these flu xes are usually far 
abo ve the fit to the SED (s ee e.g. Chini et al. i200lh . Barsony et 
al. Il998t) ). Barsony et al. il998h argue that this excess mid-IR 
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Table 3. Integrated far-infrared fluxes for all detected objects, measured with ISOPHOT, and IRAS, as well as SCUBA and IRAM 30-m points from the 
literature. Columns C50..C200 give the fluxes above the background measured with aperture photometry, columns C^q^ and C^'qq the fluxes obtained by PSF 
fitting. For comparison we list the IRAS fluxes at 12, 25, 60, and 100 /im in columns Ii2..1iO()- S450 and Sgso give SCUBA fluxes at 450 and 850 fim from 
the literature, liaoo is the flux at 1.3 mm. All fluxes are in lansky. In the B6o- B20() columns we list the background level in the ISOPHOT maps at 60, 100, 
160, and 200 fim in MIy sr^^. The * signs mark the additional detected objects in our maps. 



Object 


Ceo 


Cioo 


C16O 


C200 


"-60 


■-lOO 


Beo 


Bioo 


B16O 


B200 


S450 


Ssso 


I1300 


I12 


I25 


leo 


Iioo 


RN015FIR 


25.7 


51.5 


46.7 


48.9 


24.8 


55.5 


19 


25 


92 


120 


9.2(5) 


1.4(5) 




0.25(1) 


3.4 


47.1 


93.6 


RN015* 


3.3 


11.9 


7.1 


10.8 


3.7 


10.5 


19 


25 


92 


120 








0.18(1) 


4.2 


48.8 


82.2 


HH211-MM 


3.0 


33.8 


56.2 


54.8 


1.9 


20.8 


42 


138 


308 


331 


16.4(5) 


3.8(5) 


0.9(2) 










HH211FIRS2* 


2.9 


13.7 


55.3 


48.2 


1.0 


11.3 


42 


138 


308 


331 
















L1157 


6.8 


37.8 


42.0 


38.6 


5.9 


35.8 


11 


19 


63 


72 


6.0(6) 


0.9(6) 


0.4(6) 


0.25(1) 


0.25 


10.9 


53.5 


IC 1396W 


5.7 


19.9 


36.6 


26.6 


4.6 


13.8 


25 


92 


210 


230 








0.25(1) 


0.6 


9.7 


38.3 


IC 1396WFIRS2* 


0.1 


1.2 


33.2 


9.9 


0.2 


1.2 


25 


92 


210 


230 
















IC1396WFIRS3* 


0.1 


1.2 


4.6 


17.3 


0.2 


2.3 


25 


92 


210 


230 
















L1211 


12.9 


36.5 


63.0 


75.0 


10.5 


23.4 


42 


157 


254 


280 






0.135(3) 


2.7 


5.7 


19.8 


63.6(1) 


L1211FIRS2* 


3.4 


17.0 


31.6 


46.2 


2.2 


10.9 


42 


157 


254 


280 






0.345(4) 










CepE 


55.2 


125.8 


102.1 


81.4 


65.4 


123.0 


35 


99 


234 


293 


43.7(6) 


4.1(6) 


1.0(6) 


0.43 


5.8 


61.0 


112.0 



(I) up per lim it (2) measured at 1.1 mm (McCaughi-ean et al. 1 1994ft ') (3) mea sured a t 1.2mm (Tafall a et al ll999l) '). MMS 4 (4) measured at 1.2mm (Tafalla 
et al. Il999t) l superposition of MMS 1, MMS 2 and MMS 3 (5) Rengel et al. l200ll) (6) Chini et al. l200ll) 



37 - 



36 



H 35 - 



34 - 



33 



-1 I I I I I I I I I I r 



_L 



2.6 



2.8 



3.0 3.2 
Sn[10-i°sr] 



3.4 



3.6 



Figure 1. (T,SC) plane for the graybody fit of the CepE photometry. The 
fit with the measured fluxes is marked by a cross. Circles indicate fits using 
fluxes which deviate at most by 0.5 cr, and the small dots represent fits using 
fluxes with a maximum deviation of l.Ocr from the measurements. 



emission is due to ongoing outflow/dust interactions. Another ef- 
fect, leading to the detection of these extremely red sources at such 
short wavelength (and especially at 12 /im) could be the tiniest of 
red leaks of the IRAS filters. Suppose that there was a red leak 
of only 0.1 % in the 12 fim filter - at the limit of which the filter 
transmission curve is known - over the band pass between 20 and 
30 fjm. The IRAS detectors were still sensitive in this range. Such 
leak would increase the measured I12 flux of Cep E, for example, 
by a factor of 500. We also note that the 12 and 25 fim IRAS data 
points, if used, for most deeply embedded sources show huge devi- 
ations from the best-fitting SEDs, compared to the presumed preci- 
sion of the flux calibration of IRAS. Less red sources, like Class 1 
or 2 objects, would hardly suffer from such leaks, because their 



SEDs are a lot less steep in the 12 and 25 ^m range. Therefore, we 
decided not to use the IRAS 12 and 25 ^m data points in our SED 
fits for our very red objects. 

If one uses the measurements for the determination of the 
bolometric temperature including also the 12 and 25 /xm IRAS 
points, nevertheless, we get slightly higher values (about 4K). 
Since it is not clear how much of this small effect is due to filter 
leaks, outflow/dust interactions, or envelope emission, we do not 
list these values in Table|4] 

Our best fitting results are given in Table|4] together with the 
fit errors. In almost all cases it was not possible or not useful to do 
a fit for the newly detected objects in our maps, since they are at 
the edge of the map and so we are missing an unknown part of their 
flux. Some objects are detected only with the C200 detector since 
they are outside the slightly smaller maps at the CIOO wavelengths. 

A determination of the fit errors cannot be obtained analyt- 
ically. Therefore, we varied the measurements within their one 
sigma error box (five equidistant values; 5™ ± n/2 • AS*™; n = 
0, 1,2) and computed the best fitting parameters for each of the 
5^ = 78125 combinations. This results in an area of the parameter 
space into which the error boxes are mapped. As an example, we 
show in Fig.Qthis area in the (r,EQ) plane for CepE. The errors 
given in Tableware read off such diagrams for each of our objects. 
This procedure was applied also to the parameters L(,oi , l^smm and 

Tbol- 

The determined graybody fits are integrated to obtain the to- 
tal luminosities of the sources. By integrating only at wavelengths 
larger than 350 ^m we obtain the sub-mm luminosities hsmm, 
which can be compared to the total luminosities Lj,oi to decide 
whether an object is a Class source (Andre et al. i200(J) y Both 
values, Lfcoi and LsmmfLtoi, are given in Table|4] When the ratio 
Lsmm/Lfcoi exceeds 0.005, then the object is counted as Class 0. 
This is equivalent to th e mas s ratio Me„„/M, being larger than 
unity (see Andre et al. i2000h and references therein). The given 
bolometric temperatures Ttoi are the temperatures of a blackbody 
with the same mean frequency as the graybody, where the mean 
frequency V of an SED is determined by 
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SED{u) dv j SED(i/) dv (6) 



3 RESULTS 

Our observations of CepE, HH211-MM, IC 1396W, LI 157, 
L 1211, and RNO 15 FIR were carried out at their nominal IRAS 
positions. In our ISOPHOT maps (shown in Figs.|2|-|8) we de- 
tected more objects than were actually targeted. In four cases other 
(partly) unexpected embedded objects or bright diffuse continuum 
emission are found. For L 1211 no object was detected at the nom- 
inal IRAS position, but there were two other sources discovered 
in the maps. Measured fluxes in all filters, including IRAS fluxes 
and sub-mm and millimeter points from the literature, are given in 
Table|3l 

Discrepancies of the fluxes between PSF and "aperture" pho- 
tometry are for various reasons: First, it is a major problem to de- 
termine which pixel contributes to which object when doing "aper- 
ture" photometry. A second problem is the determination of the 
background. When using the PSF fitting method, the background is 
determined automatically (provided that the object is a point source 
and in the centre of a pixel), while for "aperture" photometry one 
has to choose background pixels. Concerning the absolute calibra- 
tion errors for the two detectors of 15 and 10% for the CI 00 and 
C200 detector, respectively, and an additional error of 20% due 
to background uncertainties, we find that both flux determination 
methods lead to consistent results in almost all cases. 

Most of our investigated objects are of Class type according 
to the LsmTn/Lfcoi criterion. We cannot decide whether the newly 
discovered objects in our maps are of Class 0, because they are 
situated at the edges of the ISOPHOT maps. Due to the different 
sizes of the maps we certainly underestimate their fluxes at 60 and 
100 /^m, which alters their derived SED in the way that they seem 
to be proportionally brighter at the longer wavelengths, but to an 
unknown extent. 

The ISOPHOT and the IRAS fluxes at 60 and 100 ^im are 
consistent within the errors only for CepE. For all other objects 
the IRAS point source catalogue gives values which are a factor of 
about 1.8 brighter. Apart from the fact that the errors for the IRAS 
data are quite large and in some cases only upper limits are given, 
the main reason for the differences is that the resolution of the IRAS 
satellite was not sufficient to resolve close-by sources. Only Cep E 
and L 1157 seem not to have other young objects in their imme- 
diate vicinity, and these are the two objects where the IRAS and 
ISOPHOT fluxes m atch th e best. Cep E is a known double source 
(Moro-Martm et al. 1200 lt» which cannot be resolved by IRAS nor 
ISOPHOT. For these reasons and the still fairly large errors in the 
flux measurements, no investigation of the time evolution of the 
fluxes of these young sources over the 14 year time span (~ 0.1% 
of the age of our objects) between IRAS and ISO is possible. 

The PSF photometry suggests that all the objects are seen as 
point sources for the ISOPHOT detectors. When subtracting the 
fitted PSF, no systematic residuals are visible in the difference im- 
ages. Thus, the angular size of the sources is at maximum 10", a 
quarter of the FHWM of the PSF. This leads to an upper limit for 
the source solid angles of about lOOD". This fact is supported by 
the inferred sizes Rioo on the order of 1 x 10"^" sr (4 □"), which is 
less than one percent of the pixel size of the ClOO detector. 

For the three objects for which we obtained a LWS spec- 
trum, we can compare the PHOT flux with the LWS continuum. 
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Figure 2. ISOPHOT maps and derived photometry for CepE. Stars (*) 
show the LWS spectrum and the solid line the best graybody fit (T = 35.3 K, 
/3= I.O, EC = 3.0 Xl0~i° sr) to the data, tioo is fixed to unity. Error bars 
for the PHOT data only include the 15% and 10% errors of the detectors. 
They do not include additional uncertainties due to background determina- 
tion. The four maps are centred on the central wavelength of the appropriate 
filter and have the same scale and orientation, given by the scale and the ar- 
row. 



While the LWS continuum is a sum of the continuum of the source 
and background radiation, the PHOT maps give the true flux of 
the source. So, the difference between LWS and PHOT should be 
the background radiation (e.g. from cold dust). In all three cases 
(Cep E, L 1 1 57, and HH 211 -MM) we clearly see evidence for such 
a background emission (see Figs.|2l-|5}. 

In the following subsections we discuss details of the results 
for the individual objects. 



3.1 CepE 

Cep E is the brightest object in our sample. Our ISOPHOT maps at 
the four wavelengths of 60, 100, 160, and 200 /im are shown in the 
lower part of Fig.|2| all at the same scale and orientation. Photome- 
try from these maps and the LWS spectrum, are displayed above the 
maps. In addition, we plot the best-fitting graybody to these data as 
solid line. The fit was used to deconvolve the measurements and the 
filter transmission curve for converting the measured fluxes to flux 
densities at the central wavelengths of the used filter. For CepE, 
the fluxes determined with PSF and "aperture" photometry were 
consistent. Deviations of the LWS continuum from the ISOPHOT 
data exist for wavelengths shorter than 100 and longer than 150 /im. 
This might be evidence for warm and cold dust. 

The PSF photometry shows that the object is a point source, 
perfectly aligned in the middle of our map, and no other embedded 
object is detected. Nevertheless Cep E is at least a double source, 
separated by 1.4" (« 1000 A U), as shown by the 222 GHz obser- 
vations of Moro-Martm et al. i200lh . 

CepE was observed by Chini et al. <200lh with SCUBA (450 
and 850 /im) and the IRAM 30-m telescope (1.3 mm). They mea- 
sured the fluxes in an aperture with a radius of 40", comparable to 
the size of our ISOPHOT pixels. The fluxes are given in Table|3l 
We included these data in the graybody fit. The deduced tempera- 
ture is 35.3 K for the best graybody fit, which also gives /3 = 1 .0 and 
ESI = 3.0x 10^^" sr. These parameters are computed fixing rioo 
to unity. If we vary the optical depth also, the rms of the fit is 
lowered from 0.9 to 0.7. The new graybody parameters are then: 
T = 42.8K, P = \.5, EQ=0.8xlO-^°sr and tioo=25. Neverthe- 
less, the inferred bolometric luminosities and bolometric tempera- 
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Table 4. Best graybody fit results, as well as the inferred bolometric and sub-mm luminosities, from our ISOPHOT data and (if available) SCUBA and 
millimeter measurements from the literature. T is the fitted temperature of the graybody. /3 the sub-mm slope of the SED, and Ef2 the solid angle of the source. 
The optical depth at 100 /im was fixed to 1.0, since it did not show significant influence on the shape of the graybody. Rioo gives the con'esponding radius of 
the envelope where tiqo is unity, determined with Tt,oi and following the assumptions of Myers et al. (1998). A discussion of fits with variable tiqo can be 
found in Sect.|3] The rms gives the deviation of the fit from the measurements scaled with the eiTors of the measurements (see Eq.|5}. The explanation of the 
determination of the errors is given in the text. Ef2 is given in 1 X 10^^*^ sr (equal to 4.25 □"). The sub-mm luminosity hsmm is the luminosity of the object 
at wavelengths larger than 350 /im, and the bolometric temperature T^oi is the temperature of a blackbody with the same mean frequency as the object. The * 
sign marks a newly discovered object. Due to the photometry problems with this object, we do not present en'ors here. 

Object T[K] 13 T.Q Rioo [AU] rms R[pc] Ttoim UodUs] Un,n,/Lboi Class 



RN015FIR 


34.0±3.0 


1.1 ±0.3 


1.7±0.2 


270 


1.0 


350 


44.6±3.0 


8.4±1.0 


0.017±0.007 


✓ 


HH211-MM 


21.0±3.0 


1.5±0.6 


12.5±1.0 


520 


2.0 


315 


31.4±1.0 


4.5±0.5 


0.046±0.020 


✓ 


LI 157 


26.5±1.5 


1.4±0.4 


3.3±0.3 


440 


1.5 


440 


37.8±1.5 


7.6±0.8 


0.025±0.015 


✓ 


IC 1396W 


30.0±2.0 


0.3 


1.1±0.2 


680 


1.3 


750 


32.6±2.0 


16.4±2.0 


0.059±0.010 


7 


L1211 


30.5±2.0 


0.0 


2.1 ±0.6 


280 


1.4 


725 


30.5±2.0 


33.1±4.0 


0.073±0.012 


7 


L1211FIRS2* 


26.8 


0.0 


1.7 


340 


1.6 


725 


26.9 


16.0 


0.100 


7 


CepE 


35.3±3.0 


1.0±0.3 


3.0±0.4 


750 


0.9 


730 


45.0±3.0 


77.9±10 


0.017±0.010 


✓ 
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Figure 3. As Fig.|5| but for L 1 157. The best graybody has the parameters 
T = 26.5 K, /3 = 1 .4, and EH = 3.3 X 10" 10 sr. noo was fixed to unity. 

tures are not affected. We obtain 77.9, 79.6 Lq and 45.0, 45.7 K for 
a fixed and free rioo, respectively. The hsmmfLtoi ratio is 0.017, a 
strong hint for the Class nature of this object, even if we observe 
the su perpos ition of two sources. Using the assumptions of Myers 
et al. il99Sh and Tboi we calculate an optical depth at 100 /im of 
9.4. From this we determine a radius of the protostellar envelope 
where rioo is unity of 750 AU. Thu s, the diameter (1500 AU) is in 
agreement with Moro-Martm et al. ^200 it) who found that the dou- 
ble system (separation of 1000 AU) is surrounded by a common 
envelope. 

The bolometric temperature of about 45 K is well below the 
value of 60 K given by Ladd and Hodapp Jl997h . They used the 
IRAS data (12, 25, 60, and 100 ^m) a nd an 800 /xm point to fit the 
bolometric temperature. Chini et al. i200lh could fit the 100 fim 
IRAS measurement and the 450, 850 and 1300 /im points. These 
data still do not cover the emission maximum of the source at about 
130 /^m. For an accurate determination of the temperature, how- 
ever, the position of the maximum of the SED is needed, which 
was observed here with ISOPHOT. 

3.2 L 1157 

Our PHOT maps of L 1 157, and the integrated photometry obtained 
from these maps, are presented in Fig.|3| This figure also shows 
our LWS spectrum of L 1157, as well as the graybody fit to the 



photometry. L 1 157 w as observed with SCUBA and IRAM 30-m 
by Chini et al. i200lh . They give two different measurements for 
the fluxes: one for the central source only (10" aperture) and one 
for the source and the whole envelope (a 55" by 30" elliptical 
aperture). The fluxes for the central source are given in TablejS] 

L 1157 is a point source with a derived temperature of about 
26.5 K, /3= 1.4, and Efi =4.3x 10"^° sr, under the assumption that 
rioo is unity. Varying the optical depth does not improve the fit. 
The bolometric temperature and luminosity are 37.8 K and 7.6 L0, 
while 2.5% of these account for Lsmm- This classifies L 1 157 as 
a Class source. According to the scheme of Myers et al. il998h 
we determine rioo ~ 13 an d Rioo = 440 AU. If we use the fluxes 
given by Chini et al. ( 2001) for the envelope, we get a much worse 
fit and it seems that we have a second cool component in the SED. 
Thus, the ISOPHOT data reflect the emission of the source itself 
and not the cold extended envelope. 

There are small deviations of the LWS continuum from the 
PHOT photometry over the whole wavelength range. This could 
be due to diffuse emission from warm and cold dust, or reflecting 
uncertainties in the calibration of ISOPHOT or LWS. PSF photom- 
etry shows that the object is not at the centre of our map, but rather 
shifted slightly to the east. This could be a hint for another source 
nearby or a slight mispointing of the telescope due to the limited 
accuracy of the IRAS coordinates. Since nothing is known in the 
literature about a second source, we attributed all the flux to L 1 157. 
PSF fitting to the C200 maps to confirm this was not possible, due 
to the unknown PSF for our half-pixel sampling in the north-south 
direction. Also the inferred size of the source (Ef2) does not sup- 
port the presence of an additional object. 



3.3 HH211-MM 

Our PHOT maps, derived photometry and a LWS spectrum of the 
HH 21 1 region are displayed in Fig.|5] HH 21 1-MM at the centre of 
our maps is the dominant source at 60 and 10 nm.l C 348 IR, prob- 
ably a heavil y embedded B-star (Strom et al. il974h , McCaughrean 
et al. il994l) ), is visible to its north-east (marked by a cross). At 
longer wavelengths, a very cold source HH 21 1 FIRS2 further north 
becomes visible and even dominant. It probably coincides with 
the o bject IC 348 MMS (marked by a circle), found by Eisloffel 
et al. 12003) to be the source of a newly detected outflow north 
of HH211. Thus, the fluxes of HH2I1FIRS2 given in Table|3lare 
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Figure 4. Best obtained fit to the SED using Eqil] Overplotted are tlie IRAS (12 and 25 /xm), ISOPHOT (60, 100, 160, and 200 ^m), SCUBA (450 and 
850 /xm), and millimeter (1100, 1300 ^tm) data points of CepE, RNO 15 FIR, L 1157, and HH211-MM (from top to bottom and from small to big circles). 
Tioo was fixed to unity. For the obtained parameters of the best fit see text or Tablel4l The models and datapoints are shifted for RNO 15 FIR, L 1157, and 
HH 2 1 1 -MM by one, two and three orders of magnitude down, respectively, for convenience. HH 211 -MM was not detected by IRAS . 
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Figure 5. As Fig.|2| but for HH211-MM. The best graybody has the pa- 
rameters T = 21.0K, /3 = 1.5, and SQ = 12.5x IQ-l" sr noo was fixed to 
unity. The position of IC 348 IR is indicated by a -I- sign, and a circle marks 
the source IC 348 MMS, found by Eisloffel et al. (2003). 



tween the PHOT photometry and the LWS continuum. So, it is very 
difficult to determine the background and to state which pixel con- 
tributes to the flux of which object. This is further complicated by 
the fact that HH211 FIRS2 has a higher surface brightness than 
our point source HH2I1-MM at 160 and 200 /im. Additionally, 
IC 348 IR could influence our measured flux for HH 21 1-MM also. 

Due to the difficulties in the determination of the fluxes of 
HH21I-MM we supplemented the ISOPHOT data with SCUBA 
datapoints at 450 and 850 ^m from Rengel et al. <2001i> and JCMT 
bolometry at 1.1 mm from McCaughrean et al. |I99J). The fluxes 
are given in Table|3| We find a fit (rms = 2.0) with the following 
object parameters: T = 21.0K, /3= 1.5, and = 12.5 x 10"^° sr. In 
this case a variable optical depth (instead of a fixed value of unity) 
does not improve the fit. The inferred bolometric temperature is 
31.4K and the bolometric luminosity is 4.5 Lq, while 4.5% ac- 
count for Lsmm- This verifies the ClassO nature of HH 21 1-MM. 
According to the low bolometric temperature we find a high rioo 
on the order of 20 and a large radius for the rioo = 1 surface of 
520 AU. 



a superposition of two different objects. The Ceo and Cioo mea- 
surements are dominated by IC 348 IR, while Cieo and C200 are 
dominated by IC 348 MMS. Therefore no further investigation of 
the SED of one of these objects was possible. 

For HH 21 1-MM, we find large differences of the fluxes at 60 
and 100 /im obtained with PSF and "aperture" photometry. These 
differences are due to the other sources influencing the background 
determination. In addition, there is a lot of diffuse background 
emission present, which can be seen in Fig.|5|as the difference be- 



3.4 IC1396W 

Our PHOT maps and the derived photometry for IC 1396W are 
shown in Fig.|6| together with the graybody fit to these data. In 
the maps, two additional sources are evident. IC 1396WFIRS2, 
to the north-east of IC1396W, peaks at 160 /im, whereas 
IC 1396 W FIRS 3, to the south-west, is remarkably red, with a flux 
ratio at C200/C160 of 2.6. These additional objects were not de- 
tected with the CIOO detector, because the maps at these wave- 
lengths are slightly smaller and the objects just fall outside. Nev- 
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Figure 6. As Fig.|2l but for IC 1396 W. We do not have LWS data for this 
source. The best graybody fit has a temperature of T = 30. OK, (3 = 03, and 
En = 1.1 X 10^^" sr Tioo was fixed to unity. Extended cool dust, or a close 
group of cold sources, are seen northeast of IC 1396 W, while a very cold 
bright source appears at 200 fim to the south-west. 
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Figure 7. As Fig.|6| but for L1211. The best fit is a blackbody with 
T = 30.5K, and EH = 2.1 X lO"^'^ sr tiqo was fixed to unity. The object 
appears north of the nominal IRAS position, and a second cool source is 
detected south-west of it. See text for the identification of these sources. 



prise, the L 1211 source was not found at its nominal IRAS posi- 
tion, but is shifted a full pixel, corresponding to 45" to the north. 
A second source L 1211 FIRS2, is found in the sou th-west. Com- 
paring our ma ps wit h the work of Tafalla et al. il99% and Anglada 
& Rodriguez J2002f) we find that our object L 1211 is identical to 
MMS 4 or VLA 5 and the object L 121 1 FIRS2 seems to be a super- 
position of MMS 3, MMS 2, and MMS 1, and VLA3 and VLA 1, 
respectively. 

The best fits to both sources are blackbodies (/3 = 0.0) and 
have temperatures of T = 30.5K and T = 26.8K for L1211 and 
L1211FIRS2, respectively. We infer a solid angle ESI of the 
sources of 2.1 and 1.7x 10~^" sr. This results in bolometric tem- 
peratures of 30.5 and 26.9 K and in bolometric luminosities of 33.1 
and 16.0 Lq, for L 1211 and L 1211 FIRS2, respectively. Even if 
we get quite high Lsmm/Ltoi ratios (0.073 and 0.100), we cannot 
firmly establish the classification of the sources as Class objects. 
This is because our ISOPHOT data do not cover the emission maxi- 
mum and also L 121 1 FIRS2 is a superposition of the emission from 
several sources. 

Tafalla et al. il999h classify L 1211 as a transitional object 
between Class and Class 1 . They use the IRAS fluxes and an ad- 
ditional observation at 1.2 mm (see Table|3}. With these data, as for 
our ISO data, the maximum of the emission could not be deter- 
mined exactly. It could only be constrained to lie between 100 and 
1200 ^m. With our ISO data, we could corroborate the assumption 
that L 121 1 is of Class 0, since the maximum of the SED is at A > 
160 /xm. If we use the 1.2 mm datapoint from Tafalla et al. ^99^ 
for our analysis, we get a poor fit. Using both, the ISOPHOT and 
the 1.2 mm point, with plausible values for (3 (1.0.. 2.0), the 1.2 mm 
flux is always overestimated by about one order of magnitude. 
The same applies for L 1211 FIRS2. Accepting values of 13 = 3.0 
or higher, we determine T = 36.3, Efi = 3.0x 10"^" sr (L 121 1) and 
T = 24.8, E!^ = 5.8xlO"^°sr (L1211FIRS2). This leads to bolo- 
metric temperatures of about 90 and 46 K and bolometric luminosi- 
ties of 77 and 19 L© for L1211 and L1211FIRS2, respectively. 
The resulting hsmm/Lboi ratios of 0.001 and 0.004 would then clas- 
sify both objects as Class 1 (L 121 1 FIRS2 is still near the transition 
phase between Class and Class 1). 



ertheless, they influence the flux measurements, especially in the 
C200 wavelengths range. We do not have observations at longer 
(sub-millimeter or millimeter) wavelengths to supplement the ISO 
data. 

For IC 1396 W we find a temperature of T = 30.0K, /? = 0.3, 
and Ef2= 1.1x10^^" sr for the best graybody fit. These parame- 
ters are obtained fixing noo to unity. A variable optical depth does 
not improve the fit. The fit is poor (rms = 1.3), however, suggesting 
errors in the photometry or a source with dust at more than one tem- 
perature. With the graybody parameters we determine a bolometric 
temperature of 32. 6K, a luminosity of 16.4 L0, and a Lsmm/Ltoi 
ratio of 0.059. Concerning the difficulties in the flux measurements 
and the not available observations at longer wavelengths, the clas- 
sification of IC 1396 W as a Class source remains questionable. 
At least sub-mm observations are needed to confirm the presence 
of a spatially extended envelope, to ensure that we do not just see a 
compact disk edge-on. 

3.5 L 1211 

Our PHOT maps of L 1211, the derived photometry, and a gray- 
body fit to these data are displayed in Fig.0 Somewhat to our sur- 



3.6 RNO 15 FIR 

The ISOPHOT maps of the RNO 15 FIR region are shown in Fig.[8l 
together with the derived photometry and the graybody fit. Vis- 
ible on our maps are RNO 15 FIR in the centre, and the source 
RNO 15 to the south-east, marked by a cross. Since this source 
is warmer than RNO 15 FIR, it is prominent at the shorter wave- 
lengths, but fades considerably relative to RNO 15 FIR towards the 
longer wavelengths (due to the larger pixel size the two objects 
also merge). From higher spatial resolution sub- mm m aps at 450 
and 850 fim taken with SCUBA (Rengel et al. j2002h ) we know 
that two other objects SMS 1 and SMS2 are present to the north and 
south of RNO 15 FIR, but are merged with it at ISOPHOT resolu- 
tion (marked by a circle in Fig.|8}. They surely influence our flux 
measurements of RNO 15 FIR. 

Since the ISOPHOT measurements show a broad and not 
well defined maximum of the SED, we su ppleme nt these data 
with SCUBA measurements of Rengel et al. i200lh to determine 
more accurate source properties. The fluxes at 450 and 850 /im 
are measured in a 45" by 45" aperture (Rengel priv. commu- 
nication) and listed in Table|3| Using these data in combina- 
tion with the ISOPHOT points, we get T = 34.0K, /3=1.1, and 
ES1= 1.7x 10^^" sr with an optical depth rioo fixed to unity. The 
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Figure 8. As Fig.|6l but for RN015FIR. The best fit has temperature 
of T = 34.0K, ^=1.1, and EH = 0.7 X 10-1° sr noo was fixed to unity. 
South-east of RNO 15 FIR, the warmer source RNO 15 is detected as well 
(marked by a cross), especially at the shorter wavelengths. The positions of 
two other weak sub-mm sources (SMS 1 - north, SMS2 - south) are indi- 
cated by a circle (Rengel at al. (2001, 2002)). 



resulting bolometric temperature and luminosity are 44.6 K and 
8.4 L0. 1.7% of the luminosity is in the sub-millimeter regime, 
classifying RNO 15 FIR as a Class source. If we vary the opti- 
cal depth at 100 /^m also, the fit is improved (the rms is lowered 
from I.O to 0.6), but th e parameters Ttoi and htoi do not change. 

Davis et al. Jl997h suggested that RNO 15 FIR might be a dou- 
ble source, due to the observed wiggling in the outflow. This might 
also be indicated by the deviation of the data points from the deter- 
mined SED (see Fig.|4j- 



4 DISCUSSION 

4.1 Mass and Age determination 

Do Class objects develop into Class 1 and Class 2 protostars? To 
answer this, we wish to determine basic parameters for the Class 
protostars, such as age, surrounding mass, present mass and final 
mass. These, however, are model dependent quantities. In Fig.|9|we 
plot the locations of our seven (including L I2I1FIRS2) sources 
on the Li,oi -Ti,oi diagram (large crosses), which is the pro tostell ar 
equivalent of a Hertzsprung-Russell diagram (Myers et al. il998h ). 
Also plotted on the diagram are the data for another 37 Class 
protostars, as listed by Andre et al. i200(]l) . Two of these sources 
possess bolometric luminosities above 1000 L©, and so fall outside 
the display. Note that Class protostars possess bolometric temper- 
atures below ~ 80 K. 

The present sample contains quite powerful and cold Class 
members. Two sources lie above the location of the other ex- 
plored sources. These are L 1211 and CepE. As we demonstrate 
below, such powerful Class sources with low bolometric temper- 
ature, can indeed be included in an evolutionary model through the 
Classes 0- I -2. The large surrounding masses observed restrict 
the type of model and these objects could go on to produce high- 
mass stars. 

The model tracks plotted represent the evolution of three pro- 
tostars which end up accumulating masses of 0.2, 1, and SM©. 
The tracks were deri ved by comb ining the Unification Scheme, as 
reviewed by Smith i200(l l2002h . with the framework for proto- 
stellar envelopes presented by Myers et al. il998h . according to 
the prescription presented below. We thus determine model ages. 
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Figure 9. The bolometric luminosity-temperature diagram for the objects 
analysed here (thick Xs), the Class data from Tablelsl (filled diamonds), 
and the Class data from the review of Andre et al. (2000) (open diamonds). 
The superimposed evolutionary tracks are discussed in Sectionlsl Protostars 
evolve from right to left. Thi'ee tracks for final masses of 0.2, 1, and 5 Mq 
are displayed. The model peak accretion rate is reached at 17,000 yr, and 
the power law fall-off is oc t"^/** with time t, on a 30.000 yr timescale. The 
vertical dotted lines on the tracks mark the model ages of 20, 30, 40, 50, 
and 75 thousand years. 



present masses of the protostellar nucleus, envelope masses and 
the final stellar masses (Table|6j. The result is that the more mas- 
sive Class protostars possess large envelopes and would become 
massive stars. According to the model described here, most of the 
envelope, however, is not accreted but dispersed, if the majority 
of protostars here are to form low-mass stars. Note that alterna- 
tive scheme s have been presented by Bontem ps et al. il996l) . Sara- 
ceno et al. jl996l) . and Andre et al. i2000l) . The model envelope 
masses are in agreement with the measurements. lust for the two 
objects where we could not determine a proper bolometric temper- 
ature (L 121 1 and L 121 1 FIRS2), there is a significant deviation. 



4.2 Outflow Luminosity vs. Source Properties 

Is the luminosity of the outflows from the Class sources corre- 
lated with the properties of the sources like their bolometric lumi- 
nosity, the temperature, or mass of their envelopes? To answer this 
question we measured the luminosi ty of the outflows from 16 of 
the Class sources in Andre et al. i200(]l) and our objects, in the 
1 -0S(1) line of molecular hydrogen at 2.122 nm. D ue to the short 
cooling times (some years. Smith and Brand il99(l) ). H2 is a good 
tracer of emission of shocked gas caused by current interactions be- 
tween outflowing material and the surrounding gas. The 1 - 0S(1) 
line of H2 is usually the brightest ro-vibrational line in a spectrum 
of shocked gas and thus most easily detected. 

I n mae netohydrodyna mic m odels of Class sources (e .g. Sh u 
et al. Jl994l> . Hirose et al. h997h . and Ouyed and Pudritz Jl997h ) 
the accretion rate onto the protostar is connected to the amount 
of material injected into the outflowing jet. This material interacts 
with the surrounding quiescent gas in shocks. Hence, the luminos- 
ity of these shocks may be connected to the mass accretion rate 
and thus to the source properties. Here, we tested for correlations 
of the l-OS(l) H2 luminosity with the source bolometric lumi- 
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Table 5. Summary of the Class sources for which a correlation of the 
source properties with the outflow luminosity in the 1 - S( 1 ) line of H2 
was investigated. Except for our objects observed with ISOPHOT, L;,;,; and 
Tf,oj are adapted from Andre et al. (2000), as well as all the Menjj values 
(except RNO 15 FIR and L 1211). The outflow luminosities are either from 
published literature or our own measurements. The typical errors of our 
outflow luminosities are 10%. In the Ref. column the references are given 
where we took Lhj i-osm measurements from. 



Object 




T(,o! 




l-OS(l) 


Ref. 




[L©] 


[K] 


[Mq] 


[10-3 Lq] 






RNO 15 FIR 


8.4 


44.6 


0.9* 


0.46 


1, 


7 


HH211-MM 


4.5 


31.4 


1.5 


3.1 


2, 


7 


L1157 


7.6 


37.8 


0.5 


6.1 


3, 


7 


IC 1396W 


16.4 


32.6 




19.0 


6 




L1211 


33.1 


30.5 


0.8** 


10.7 


7 




CepE 


77.9 


45.0 


7.0 


70.0 


4, 


7 


L1448N 


11.0 


70.0 


2.3 


2.16 


7 




L 1448 IRS2 


6.0 


70? 


0.9 


2.8 


7 




L1448C 


9.0 


60.0 


1.4 


5.7 


7 




IRAS 03282 


1.5 


35.0 


0.6 


4.46 


7 




HH212MM 


14.0 


70? 


1.2 


5.3 


5, 


7 


HH24 MMS 


5.0 


20? 


4.0 


1.21 


1, 


7 


HH25MMS 


6.0 


34.0 


0.5 


6.61 


1 




NGC2264GVLA2 


12.0 


25.0 


2.0 


7.75 


3 




VLA 1623 


1.0 


<35 


0.7 


0.81 


3, 


7 


Ser-FIRSl 


46.0 


51.0 


3.0 


0.64 


7 





* taken from Rengel et al. 120031) 
** taken from Tafalla et al. < 1999ft 

References: mOavis et al. 1 19971) (2) McCaughrean et al. <1994l) (3) D avis 
& Eisloffel 1 19951) ( 4) Eisl offel et al. Il99el) (5) Zinnecker et al. Il998t) (6) 
Froebrich & Scholz J2003I) (7) own measurements 



nosity, the bolometric temperature and the m ass of the surrounding 
protostellar envelope given by Andre et al. i200Cl) . The results of 
these comparisons are shown in Figs. llOlfm and ll2l We obtained 
a linear regression for each case and tested if the slope of the re- 
gression line differed statistically significantly from a slope value 
of zero. With a probability of error of 5 % none of the regression 
lines differs from a constant value. Additionally a Kolmogorow- 
Smirnow-Test shows, that with a probability of error of 0. 1 % the 
data is not consistent with a constant value. Thus, a significant cor- 
relation of the outflow luminosity in the 1 - 0S(1) line of H2 with 
any source parameter was not found. 

The lack of such correlations may have various explanations. 
For example, in each outflow, we observe H2 emission at various 
distances from the source and these knots or bow shocks are indi- 
cating material which was ejected from the source at diiferent times 
in the past (At = distance to the source/jet velocity). Also, the knot 
luminosity depends on the local properties of the surrounding gas 
(e.g. gas density, atomic fraction). Additionally, the extinction gra- 
dient in the K-band along the outflow is not known. It will alter the 
measured relative and total fluxes in the sense that knots closer to 
the source appear fainter due to higher extinction. In Section |5] 
however, we argue that the location on these diagrams depends 
sensitively on both mass and age, which results in a wide scatter. 
Therefore, the lack of a significant correlation of the present source 
properties and the outflow luminosity in the H2 l-OS(l) line may 
not be surprising. 

A better tool for comparing the outflows to source proper- 
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Figure 10. Measured outflow luminosity in the 1 -0S(1) line of H2 ver- 
sus the bolometric source luminosity for the sources listed in Tablelsl The 
objects investigated in this paper are marked with a large + sign. No signif- 
icant congelation is found for this sample. 
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Figure 11. As Fig. llOl but for the bolometric temperatures of the Class 
sources from Tablelsl No significant correlation is found for this sample. 



ties may be an optically thin line of CO (e.g. the 1-0 CO line), 
which should give a measurement of the time-integrated power 
of the outflow without being influenced by local extinction ef- 
fects. Comparable observations in the same transition and isotope 
of CO are needed for a statistically reasonable sample of objects 
to study their behaviour. At present, only small samp les of Class 
sources have been thus analysed (e.g. Bontemps et al. l ll99dh . Smith 
i200(il2002h ). 
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Figure 12. As Fig. llOl but for the envelope masses of the Class sources 
from Tablets! For IC 1 396 W we do not have measurements of the envelope 
mass. No significant correlation is found for this sample. 



5 AN EVOLUTIONARY SCHEME 

An evolutionary model for protostars is presented in t he Ap 
The o utflow scheme has been elab orated by Smith {^98 
l2002h and a pplied by Davis et al. il99A . Yu et al. hOlWh . and 
Stanke et al. j200Cl) . It is based on a prescribed accretion rate from 
an envelope. Modelling of outflows has demonstrated that the frac- 
tion of mass which escapes through jets must reach a maximum 
during the Class stage. This is required to account for the excess 
momentum and power of Class bipolar outflows, as calculated 
from obser vations of emission lines of CO rotational transitions 
(see Smith ilOOOh ). We outline in the Appendix the fundamental 
formulae of the evolutionary scheme. 

According to previous modelling of the envelope, three pa- 
rameters must be introduced to generate plausibl e mod els for the 
bolometric temperature. As shown by Myers et al. Il998h . these are 
(1) the envelope's outer temperature (here To =24K), (2) the effi- 
ciency of accretion of the envelope into the star-jet system and (3) 
the difference in evolutionary timescale between the envelope and 
the protostar. The envelope consists of material which will fall onto 
the central object as well as mass directly lost soon after the Class 
stage. This extra mass component proves necessary to produce a 
low bolometric temperature, as observed for the Class sources, 
yet must be rapidly lost in order to yield T Tauri stars within a rea- 
sonable time (Myers et al. il998l) '). 

Previously, we modelled the envelope evolution by assuming 
mass conservation. Here, we find two significant adjustments are 
necessary in order to model the new data and maintain plausible 
time scales. First, the initial mass in the envelope which will even- 
tually fall inwards is reduced to 87 % of the total required mass 
to form the star and excavate the bipolar outflow. The other 13 % 
is presumed to initially lie within a flattened disk. This yields the 
values in the column 'infall mass' in Table|6l 

In Fig.|9| we plot the sample summarised in Table|5|for Class 
sources for which a correlation of the source properties with the 
outflow luminosity in the l-OS(l) line of H2 has been investi- 
gated. Note that this sample includes warmer and less luminous 




0.0 



0.5 



1.0 1. 
log(L,o,/Ls) 



2.0 



2.5 



Figure 13. The derived outflow shock luminosity versus the bolometric 
source luminosity for the Class sources listed in Tablelsl as well as the 
Class sample investigated by Stanke (2000) (symbol: 'x') and the far- 
infrared line ISO luminosities presented by Giannini et al. (2001) (symbol: 
'*'). The model tracks are for the same three models presented in Figl9l 
and the straight line divides model Class and model Class 1 protostars, as 
determined by the protostar possessing half of its final mass. 
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Figure 14. The derived outflow shock luminosity versus the bolometric 
source luminosity for Class 1 sources from the sample of Stanke (2000) 
(symbol: 'x') and the far-infrared line ISO luminosities presented by Gi- 
annini et al. (2001) (symbol: '*'). The model tracks are for the same three 
models presented in Fi gl9l and the straight line divides model Class and 
model Class 1 protostars, as determined by the protostar possessing half of 
its final mass. 



protostars than in the ISOPHOT sample investigated above. Ac- 
cording to the tracks, this corresponds to a wide range in final stellar 
masses. The lowest mass star forming here is found to be VLA 1623 
(assuming Tjoj = 35 K) which will reach just 0.07 M© , owing to its 
low bolometric luminosity of only I Lq (Andre et al. |2000|)). The 
low final mass is a result of this version of the evolutionary scheme 
employed, for which we maintain the same accretion timescale but 
alter the accretion rate to generate the tracks. This implies that the 
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Table 6. Parameters derived for the seven objects from the model evo- 
lutions. The minimum mass is the total mass with density distributed as 
p oc r"-^/^ necessary to provide an optically thick sphere out to a radius 
Rjjoj, coiTesponding to the observed Ttol- The infall mass is the envelope 
mass which remains to be accreted (a part of which will escape in the jets), 
and the envelope mass is the total mass predicted on projecting the dis- 
tribution out to a radius corresponding to the chosen ambient temperature 
of 24 K. The model mass accretion rate decreases as t~ on a 30,000 yr 
timescale. The age is given in 10'' years, the masses are in solar masses. 
For comparison we give the measured values for the envelope masses from 
the literature in column Menv- The values correspond to the Menv column 
in Tablelsl The envelope mass for L 121 1 FIRS2 is taken from Tafalla et al. 
(1999). 



Object 


Age 


Mass 


Final 
mass 


Min. 
mass 


Infall 

mass 


Env. Menv 

mass 


RN015FIR 


29.7 


0.10 


0.5 


0.05 


0.32 


1.0 


0.9 


HH211-MM 


24.0 


0.06 


0.3 


0.19 


0.26 


1.4 


1.5 


L1157 


27.0 


0.10 


0.5 


0.10 


0.35 


1.3 


0.5 


IC1396W 


24.9 


0.21 


1.2 


0.50 


0.88 


4.2 




L1211 


24.0 


0.43 


2.6 


1.53 


1.91 


10.3 


0.8 


L1211FIRS2 


22.5 


0.21 


1.4 


1.42 


1.06 


6.9 


2.1 


CepE 


30.6 


0.97 


4.2 


0.37 


2.83 


8.1 


7.0 



final mass is nearly proportional to the peak accretion luminosity. 
Future statistical studies will lead to revisions of this first model. 

The simplest form of the unifying model, assumed here, is that 
a fraction of the jet power is instantaneously dissipated in shock 
waves, while the bipolar outflow is a time-averaged recording of the 
momentum outflow. To model the outflow, we have previously em- 
ployed the H2 luminosity, L(H2), which we estimate to be ten times 
the 1 - 0S(1) luminosity. This is consistent with expectations from 
shock physics and allows a comparison with previous diagrams pre- 
sented by Stanke ( 2000) and Smith Here, however, we shall 
use the jet power itself as the comparison parameter. For the com- 
parison, we assume that the observed emission is produced in the 
warm shocks where the jets dissipate their energy, hghock- 

The fraction of the jet power dissipated in molecular hydrogen 
lines is taken to be 2 %. This is consistent with numerical simula- 
tions and bow shock modelling which predict, typically, 10% of 
the infrared radiation from shocks in dense clouds to be in the form 
of H2 lines. We also assume that 80 % of the jet energy is hidden by 
just under two magnitudes of K-band extinction. The shock power 
has also been estimated from the far-infrare d lines of CO, 01, OH, 
and H2O, measured by ISO (Giannini et al. J200lh ). Here, we shall 
assume that these lines in total, within the ISO-LWS beam, also 
represent 2 % of the jet power, hshock- We thus increase estimated 
H2 and sub-mm luminosities by 50 to yield the displayed values. 
While these approximations are far from ideal, they permit us to 
determine if the evolutionary scheme is plausible. 

Figure fTsl demonstrates that the Class protostars possess al- 
most exclusively high ratios of 'L(shock)fLboi- Two objects, how- 
ever, RNO 15 FIR and Ser-FIRS 1, lie well within the model Class 1 
regime. Environmental factors could cause a downward shift of the 
data points, including higher extinction or radiative shocks which 
are less efficient in H2 vibrational excitation. Furthermore, these 
data are consistent with the same model tracks fitted to the bolo- 
metric luminosity-temperature data. The lack of a correlation in 
the data is thus put down to the combination of the distributions in 
both mass and age. 

In addition, previously measured Class 1 outflows almost all 
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Figure 15. The derived outflow shock luminosity versus the bolometric 
source temperature for the Class sources listed in Tablelsl The model 
tracks ai'e for the same three models presented in Fi °l9l 

lie below the predicted Class border line, as shown in Fig.^| Note, 
howev er, that for many of the H2 flows in Orion detected by Stanke 
("2000') only upper limits for the bolometric luminosity are avail- 
able. Nevertheless, the division of the two Classes with the model 
straight line where the protostar has acquired half the final stel- 
lar mass, is evident. A similar difference in outflow luminosity be- 
tw een Cl ass and Class 1 sources has been found by Bontemps et 
al. <199d) . 

Figures fTsl demonstrates that the model is also consistent with 
the envelope properties. The main exceptions apparent from this 
diagram are a group of low luminosity H2 objects. This suggests 
that the extinction for these sources may far exceed the fiducial two 
magnitudes. 



6 CONCLUSIONS 

We have observed the spectral energy distributions for the six 
deeply embedded objects CepE, L1211, IC 1396W, L1157, 
HH211-MM, and RNO 15 FIR in the far-infrared with ISO. The 
inferred temperatures and hsmmfl^boi ratios confirm the Class na- 
ture of four of these sources, within the errors. Employing an evolu- 
tionary scheme, we are able to estimate the age, surrounding mass 
and the current and final mass of these sources. These estimates are, 
however, model dependent. Two sources, CepE and L 121 1 appear 
to develop into intermediate mass stars, while the others will be- 
come solar mass stars or lower mass objects. The comparison of 
the ISOPHOT and LWS observations for three of the sources re- 
veals the existense of emission from cold dust in the immediate 
vicinity of the objects. 

A comparison of the luminosity in the 1 - 0S(1) line of H2 of 
the related outflows for 16 Class sources, with the source bolo- 
metric luminosity, bolometric temperature, and envelope mass was 
done. We found no statistically significant correlation of the out- 
flow luminosity with each of these source parameters. This could 
be due to the H2 luminosity mainly depending on the local proper- 
ties of the surrounding gas. 

The unifying scheme, however, explains the lack of correla- 
tions as due to evolutionary effects. Furthermore, the scheme which 
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involves a redistribution of mass between envelope, disk, protostar, 
jets and outflow, accounts for the differences in source properties 
according to the Class. 
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APPENDIX 

We present and test a model based on the transfer of gas between 
components. We take a spherical envelope of gas and dust, and pre- 
scribe an accretion rate from the inner edge of the envelope onto a 
disk. Note that we assume a centrifugal bamer at 30 AU, which de- 
fines the inner envelope - outer disc transition. The accretion disk 
processes most of the mass onto the protostar and a fraction into 
twin jets. The speed of the jets is assumed to be a fixed fraction of 
the escape speed from the protostellar surface. 

The accretion rate from the envelope is taken to increase ex- 
ponentially for a short period before decreasing as a power law 
through the Class 0, 1 and 2 phases. The zero point of time is thus 
defined as the moment when accretion starts and, simultaneously, a 
central hydrostatic object forms. The accretion rate is 



Ma{t) = Mo{e/a)°'{t/to)-°' exp(-to/t). 



(7) 



Energy release through accretion and contraction are included. In 
the models shown, the peak accretion rate is reached at to/a = 
17, 000 yr, and the power law index is a = 7/4, on a to = 30,000 yr 
timescale. The accreted mass is predominantly accrued by the 
growing protostars. The fraction e{t) which escapes through twin 
jets reaches a maximum of 77 = 0.2 at the peak accretion time: 



Ma[t) 
Mo 



(8) 



where = 2 is found to be appropriate. Hence the mass left over, 
which accretes onto the core to form the star is 



(1 - €)Ma. 



(9) 



To form a star like the Sun, this model will provide an early 
accretion peak in which Ma ~ 1O~*M0 yr~^ for 10^ years, and 
eventually becoming Ma ~ 1O~^M0 yr~^ for 10*^ years, corre- 
sponding to Class and Class 2 or Classical TTauri stars, respec- 
tively. Thepower-law has substantial observational support (Calvet 
et al. ilOm ). 

We previously modelled the envelope evolution by assuming 



mass conservation. Here, we make two significant adjustments in 
order to model the new data. First, the initial mass in the envelope 
which will eventually fall inwards is reduced to 87% of the total 
required mass to form the star and excavate the bipolar outflow. 
The other 13 % is presumed to initially lie within a flattened disk. 
This yields the values in the column 'infall mass' in Table|6| The 
total mass can be written analytically in terms of an incomplete 
Gamma function on integrating Eqn.0 



M,„faii{t) = Moto{e/a)°' [1 ~ r{a - l,ta/t)] . 



(10) 



Secondly, we find that the low bolometric temperatures of Class 
protostars can only be attained by introducing an additional mass 
compone nt to th e envelope. In confirmation of the results of My- 
ers et al. il998l) . we find that this extra mass is lost on a shorter 
timescale than the protostellar accretion timescale. The bolometric 
temperature is calc ulated according to the optically thick case of 
Myers et al. il998h . We thus extend the opacity law approximation 
employed up to 60 to 120 fim with the same form and take the op- 
tically thick envelope throughout the early evolutionary stages. We 
have thus found here that an envelope mass 



Me„4f) = M,„f ■ (0.87 + /i(Vio)"'") 



(11) 



where /j. = 2 provides bolometric temperatures, timescales and 
masses consistent with the observed samples. 

The envelope mass provides a testable prediction. This mass 
is not strongly dependent on the evolutionary path but is neces- 
sary to provide the optical depth out to a sufficiently large radius 
to permit the measured low bolometric temperature. The total mass 
is dominated by the outer regions of the envelope, while the to- 
tal optical depth is controlled by the inner region (for all plausible 
density distributions such as p oc r"'^''^, as assumed here). Hence, 
the mass is sensitive to the extent of the envelope. For this rea- 
son, we present three determinations of the envelope mass in Ta- 
ble|6| Masses derived from submill imetre observations yield quite 
low extended masses (Andre et al. i2000h ). consistent with the ab- 
sence of m ore ma ss than necessary to form the star and feed the 
jets (Smith i200Cl) ). It is clear that both the observationally derived 
mass and model mass are sensitive to chosen physical parameters 
and both will need refining. 
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